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Abstract 


Predicted  performance  for  antenna-coupled,  quantum-well  infrared 
photodetectors  (QWIPs)  is  presented.  These  structures  contain  a  pat¬ 
terned  metal  antenna  layer  and  dielectric  backplane  reflector  in  place  of 
the  usual  metal  phase  grating.  When  illuminated  by  normally  incident 
plane-wave  illumination,  the  antenna  structures  generate  both  evanescent 
and  propagating  modes  that  have  the  required  polarization  to  be  ab¬ 
sorbed  by  the  QW  stack.  An  undersampled  spectral  technique  is  also 
demonstrated  for  using  the  finite-difference,  time-domain  (FDTD)  tech¬ 
nique  to  compute  spectral  quantum  efficiency  in  a  single  FDTD  run. 
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Introduction 


Quantum-well  infrared  photodetectors  (QWIPs)  provide  many  attractive 
features  for  use  in  large  focal  plane  arrays  (FPAs).  Most  QWIPs  are  grown 
on  GaAs  substrates,  for  which  large  wafers  (up  to  6  in.  diameter)  are 
available.  Reproducible  and  precise  growth  using  established  III-V  mo¬ 
lecular  beam  epitaxy  (MBE)  and  a  simple  fabrication  process  has  pro¬ 
duced  high  yields  in  both  detectors  and  FPAs.  The  large  wafers  and  uni¬ 
formity  that  can  be  achieved  with  MBE  offer  a  unique  capability  for  low- 
cost  production  of  large,  high-sensitivity  staring  FPAs.  640  x  480  FPAs 
with  thermal  sensitivity  around  20  mK  are  now  typical. 

However,  the  quantum  efficiency  (QE)  of  QWIP  detectors  is  typically 
lower  than  that  of  direct  bandgap  detectors,  such  as  HgCdTe  photo¬ 
diodes.  A  typical  QWIP  QE  is  in  the  range  of  10  to  30  percent,  while  direct 
bandgap  detectors  achieve  70  to  80  percent.  The  spectral  bandwidth  of 
QWIPs  is  also  typically  limited  to  -1.5  pm,  further  reducing  the  detection 
efficiency  for  blackbody  targets.  For  many  applications  in  which  the  de¬ 
tection  sensitivity  is  limited  by  the  charge  well  capacity  of  the  readout  cir¬ 
cuitry,  the  difference  in  QE  has  no  effect.  However,  in  other  applications — 
e.g.,  systems  operating  at  high  frame  rates  or  in  low  backgrounds — 
higher  QE  and  bandwidth  would  be  beneficial. 

A  typical  QWIP  pixel  is  shown  in  figure  1.  It  consists  of  a  QW  stack  sur¬ 
rounded  by  two  contact  layers  and  an  optical  coupling  structure  (OCS), 
such  as  a  metalized  diffraction  grating.  The  OCS  is  required  to  scatter  or 
diffract  normally  incident  light  into  modes  with  a  component  of  electric 
field  perpendicular  to  the  QWs.  The  efficiency  of  this  scattering  or  diffrac¬ 
tion  process,  and  the  degree  to  which  the  scattered  light  is  trapped  inside 
the  pixel  to  make  multiple  passes  through  the  QW  stack,  are  major  factors 
that  determine  the  QE  of  the  QWIP 


Figure  1.  Typical 
grating-coupled  QWIP 
pixel. 
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Previous  QWIPs  have  used  a  variety  of  metalized  diffraction  gratings,  in¬ 
cluding  periodic  [1-5]  and  pseudorandom  [6]  etched  gratings,  as  well  as 
periodic  planar  gratings  [7],  These  structures  are  generally  designed  to 
couple  via  propagating  modes  above  the  grating.  With  periodic  gratings, 
the  period  is  chosen  so  that  the  grating  is  near  resonance,  i.e.,  so  that  nor¬ 
mally  incident  light  is  diffracted  at  an  angle  slightly  less  than  90°  from  the 
normal.  In  a  large  pixel,  the  diffracted  light  makes  a  pass  through  the  QW 
stack,  is  totally  internally  reflected  by  the  GaAs/air  boundary,  makes  a 
second  pass  through  the  QW  stack,  and  is  then  coherently  rediffracted 
back  out  of  the  sample  when  it  reimpinges  on  the  grating  [6].  This  coher¬ 
ent  rediffraction  happens  when  all  of  the  diffracted  orders  have  the  same 
path  length  between  leaving  and  returning  to  the  grating,  so  that  they  re¬ 
turn  with  the  same  phase  relationship  with  which  they  left.  The  pseudo¬ 
random  grating  scatters  normally  incident  light  into  a  wider  range  of 
angles  so  that  the  totally  internally  reflected  modes  reimpinge  on  the  grat¬ 
ing  with  randomized  phase,  thereby  avoiding  coherent  rediffraction  out 
of  the  QWIP,  and  permitting  additional  passes  through  the  QW  stack.  In 
large  pixels,  a  significant  improvement  of  QE  was  observed  with  the 
pseudorandom  over  the  periodic  grating.  However,  the  improvement  de¬ 
creases  in  smaller  pixels  [8].  The  reduction  in  QE  for  smaller  pixels  tends 
to  be  a  characteristic  of  all  QWIPs  that  rely  on  propagating  modes,  and  is 
a  problem  for  the  small  (typically  ~24  pm)  pixels  used  in  large  arrays. 
Periodic  gratings  also  couple  effectively  in  a  relatively  narrow  (1-2  pm) 
spectral  band.  Pseudorandom  gratings  typically  have  a  wider  spectral 
band  corresponding  to  the  wider  range  of  spatial  frequencies  in  their 
patterns. 

Several  groups  have  demonstrated  designs  in  which  the  "grating"  is 
etched  through  the  QW  stack,  either  in  a  two-dimensional  square  lattice 
(the  "enhanced,"  or  EQWIP)  [9]  or  as  a  set  of  parallel  V-shaped  grooves 
(the  "corrugated,"  or  C-QWIP)  [10].  These  designs  depend  less  on  propa¬ 
gating  modes  and  more  on  localized  modes,  so  their  performance  is  de¬ 
graded  much  less  in  small  pixels.  The  V-groove  structure  [10]  also  has  a 
relatively  wide  spectral  bandwidth. 

In  this  report,  we  investigate  the  expected  performance  of  antenna- 
coupled  QWIPs.  These  structures  employ  an  array  of  planar  metal  anten¬ 
nas  in  close  proximity  to  the  QW  stack,  so  that  the  near-field  modes  of  the 
antennas  couple  to  the  QWs. 

Design  and  modeling  of  optical  coupling  in  QWIPs  is  difficult  because  the 
absorption  takes  place  in  the  near-field  of  the  grating  (invalidating  the 
common  Presnel  and  Praunhofer  approximations),  the  features  are  com¬ 
parable  to  the  wavelength  of  the  optical  radiation,  and  the  anisotropic  ab¬ 
sorption  of  the  QW  stack  requires  use  of  a  vector  propagation  model  (ex¬ 
cept  in  certain  simple  geometries).  Eor  structures  in  which  the  electro¬ 
magnetic  modes  in  the  grating  are  known,  such  as  the  linear  lamellar 
grating  or  the  two-dimensional  array  of  box  cavities,  the  coupling  can  be 
calculated  by  the  modal  expansion  method  (MEM)  [4].  However,  the 
modes  are  known  for  only  a  few  simple  structures  and,  even  then,  only 
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for  the  case  of  perfect  electric  conductor  (PEC)  grating  metal.  For  ex¬ 
ample,  while  it  is  possible  to  model  a  two-dimensional  array  of  square 
cavities,  the  complementary  structure  with  square  pedestals  is  much 
harder. 

Previous  experimental  results  and  the  results  presented  later  in  this  report 
show  that  finite  conductivity  of  the  metal  used  for  the  grating  or  antenna 
can  significantly  degrade  QE.  For  example,  earlier  experiments  [4]  indi¬ 
cated  that  detectors  with  AuGe  contacts  had  a  QE  that  was  70  percent  of 
the  QE  for  detectors  with  pure  gold  contacts.  The  results  in  this  report 
suggest  that  the  finite  conductivity  of  even  pure  gold  is  significant. 

This  report  describes  the  use  of  the  finite-difference,  time-domain  (FDTD) 
method,  as  well  as  a  modified  MEM  for  calculating  optical  QE  in  QWIPs. 
An  undersampled  spectral  FDTD  technique  was  used  that  is  especially 
advantageous  for  QWIP  modeling  because  it  permits  calculation  of  a  full 
spectrum  of  QE  versus  wavelength  in  a  single  FDTD  run.  The  three- 
dimensional  absorption  pattern  can  also  be  visualized  at  any  of  the  wave¬ 
lengths  in  the  spectrum.  The  modified  MEM  is  used  to  model  QE  in  infi¬ 
nitely  periodic,  antenna-coupled  QWIPs. 

Antenna-Coupled  QWIP  Structure 

The  basic  elements  of  an  antenna-coupled  QWIP  are  shown  in  figure  2.  In 
addition  to  the  standard  QW  stack  with  common  and  top  contact  layers 
and  antireflection  coating,  the  structure  contains  a  planar  metal  antenna 
layer,  backplane  dielectric,  and  backplane  metal.  In  antenna  terminology, 
the  backplane  dielectric  and  metal  serve  to  reflect  the  back  radiation 
modes  of  the  antenna  to  the  front.  Alternatively,  they  can  be  thought  to  re¬ 
flect  light  that  passes  through  the  antenna  back  into  the  absorbing  part  of 
the  detector.  The  backplane  structure  permits  the  antenna  to  be  used  in  a 
"reflection"  configuration  rather  than  the  arrangement  used  previously 
for  transmission  planar  gratings  [7]. 


Figure  2.  Cross  section 
of  antenna-coupled 
pixel. 
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The  antenna  dimensions  are  largely  determined  by  the  wavelength 
band(s)  that  are  to  be  detected.  For  example,  the  spiral  antenna  (to  be  de¬ 
scribed  later)  fits  in  a  square  that  is  about  6  pm  on  a  side.  To  provide  cou¬ 
pling  for  the  whole  pixel,  the  antenna  pattern  is  replicated  in  a  periodic 
array  that  nearly  fills  the  pixel.  The  resulting  periodic  antenna  pattern  can 
couple  through  both  propagating  Floquet  modes  generated  by  the  peri¬ 
odic  array  and  bound  antenna  modes.  If  the  propagating  modes  make  a 
dominant  contribution  to  QW  absorption,  then  the  structure  is  similar  in 
operation  to  the  planar  metal  grating.  However,  when  the  QW  coupling 
arises  mainly  through  bound  antenna  modes,  a  potentially  broadband 
coupling  is  achieved  that  is  scalable  to  small  pixels. 

Use  of  planar  metal  layers  eliminates  the  need  for  a  grating  etch  and  the 
potential  nonuniformity  associated  with  that  etch.  On  the  other  hand,  the 
antenna  structure  requires  deposition  of  the  backplane  dielectric  and  of¬ 
ten  involves  definition  of  fine  features  in  the  antenna.  All  the  processing, 
except  for  deposition  of  the  antireflection  (AR)  coating,  is  done  on  the 
front  (MBE  growth)  side  of  the  wafer,  so  there  is  no  need  for  backside 
processing  of  individual  die  or  handling  of  fully  thinned  wafers,  as  with 
some  other  approaches. 

Modeling  Techniques 

Two  models  were  used  to  calculate  the  absorption  QE  of  QWIP  pixels. 

The  first  method  is  an  MEM,  modified  to  include  an  arbitrarily  patterned, 
finite  or  perfectly  conducting  antenna  layer.  The  procedure  assumes  the 
antenna  pattern  is  periodic  and  infinite  in  lateral  extent. 

The  second  model  is  based  on  the  three-dimensional  EDTD  method.  It 
can  model  arbitrarily  shaped,  grating-  or  antenna-coupled,  finite-sized 
QWIP  pixels.  Pinite-conductivity  metals  can  be  used,  although  with  some 
limitations. 

Coupled  Mode  (Infinite  Pixel) 

The  variation  of  the  MEM  that  was  used  for  infinitely  periodic,  antenna- 
coupled  QWIPs  was  implemented  using  the  4x4  matrix  method  de¬ 
scribed  by  Yeh  [11].  This  transfer-matrix  method  rigorously  computes 
vector  electromagnetic  propagation  of  a  plane  wave  through  a  stack  of 
uniaxially  anisotropic  planar  layers.  In  the  case  of  a  QWIP  pixel  contain¬ 
ing  a  periodic  array  of  antennas,  the  4x4  matrix  method  was  applied  to 
each  of  the  Ploquet  modes  independently,  except  at  the  antenna  layer, 
where  the  Ploquet  modes  were  mixed  to  achieve  the  required  boundary 
conditions  along  the  antenna  surface.  The  antenna  layer  was  assumed  to 
be  a  patterned  sheet  conductor  with  zero  thickness  and  a  specified  (some¬ 
times  infinite)  sheet  conductance.  The  calculation  enforced  the  condition 
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that  the  tangential  electric  fields  were  zero  along  the  metal  parts  of  the 
antenna  (for  a  perfectly  conducting  antenna  metal)  or  that  sheet  current 
(and  therefore  the  discontinuity  in  tangential  magnetic  field  across  the 
antenna  layer)  was  proportional  to  the  tangential  electric  field  through 
Ohm's  law: 


/metal“  ^metal^ll  •  (1) 

The  antenna  was  laid  out  on  a  square  grid.  For  example,  a  unit  cell  of  the 
spiral  antenna  (see  fig.  3)  was  laid  out  on  a  19  x  19  mesh.  The  19  x  19  real- 
space  mesh  was  then  used  in  a  solution  with  a  set  of  19  x  19  Floquet 
modes. 

Finite-Difference  Time  Domain  (Finite  Pixel) 

A  QWIP  pixel  is  a  somewhat  challenging  structure  in  which  to  compute 
electromagnetic  coupling.  The  structure  is  inhomogeneous,  including 
multiple  layers  of  dielectrics  and  finite-conductivity  metals.  The  QW 
stack  itself  is  anisotropic.  Finally,  the  structure  has  many  potential  reso¬ 
nances  that  result  from  the  periodicity  of  the  grating,  as  well  as  from  the 
pixel  cavity  itself. 

The  FDTD  method  [12,13]  is  well  suited  to  solving  problems  of  this  type. 
The  method  is  a  direct  solution  for  Maxwell's  time-dependent  curl  equa¬ 
tions  based  on  volumetric  sampling  of  the  unknown  field  distribution 
within  and  near  the  pixel  over  a  period  of  time.  The  sampling  in  space  is 
typically  at  10  to  20  samples  per  wavelength  (in  the  material).  The  sam¬ 
pling  in  time  is  selected  to  ensure  numerical  stability.  For  the  results  pre¬ 
sented  here,  the  QWIP  pixel  structures  were  first  designed  and  gridded 
using  a  custom  Mathematica  notebook.  The  FDTD  calculations  were  then 
run  using  a  commercial  FDTD  program  called  XFDTD  [14].  Finally,  the  re¬ 
sults  were  analyzed  using  custom  software  in  C-i-i-. 


Figure  3.  3  x  3  array  of  — 
four-arm  spiral 
antennas.  Metal  is 
dark. 
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High,  but  finite-conductivity  metals  present  a  problem  in  FDTD  (and 
other  numerical  electromagnetic  techniques).  This  is  because  the  fields 
change  so  rapidly  near  the  surface  of  such  a  metal,  thereby  requiring  an 
unreasonably  small  grid  spacing  to  accurately  describe  the  fields.  For  ex¬ 
ample,  a  grid  spacing  of  0.3  pm  is  acceptable  in  most  sections  of  a  QWIP 
being  illuminated  with  incident  radiation  at  9.0  pm  wavelength,  since 
0.3  pm  is  about  1/10  the  wavelength  in  GaAs.  However,  a  high- 
conductivity  metal  such  as  gold  may  have  a  skin  depth  of  -100  A,  requir¬ 
ing  a  grid  spacing  of  -10  A,  which  requires  far  too  many  grid  cells  to  cal¬ 
culate  in  any  reasonable  time.  Instead  of  this  "brute-force"  approach,  we 
used  the  "synthetic  conductivity"  method  [15],  which  employs  a  derived 
conductivity  value  that,  when  used  in  the  FDTD  difference  equations, 
yields  the  correct  surface  impedance.  Therefore,  the  fields  outside  the 
metal  are  nearly  correct,  even  though  the  fields  inside  may  not  be.  It 
should  be  noted  that  although  the  synthetic  conductivity  method  has 
been  shown  to  yield  the  correct  reflectance  for  planar  surfaces  [15],  its  ac¬ 
curacy  on  highly  structured  surfaces  such  as  gratings  has  not  been  proven 
yet.  Also,  accuracy  may  be  worse  for  thin  conductors  such  as  the  antenna 
that  are  less  than  five  cells  thick. 


The  appropriate  synthetic  FDTD  conductivity  is  related  to  the  true  con¬ 
ductivity  through  [15] 


^FDTD- 


(2) 


where  Az  is  the  cell  size,  co  is  the  angular  frequency  of  the  illumination, 
and  p  is  the  magnetic  permeability.  For  a  good  conductor,  we  use  the  ap¬ 
proximation  that  the  reflectivity  is  approximately  related  to  the  conduc¬ 
tivity  through  [16] 


1-R 


V  f^Tme 


(3) 


where  is  the  permittivity  of  free  space.  Eliminating  from  the 
equations  and  assuming  nonmagnetic  materials  yields  the  frequency- 
independent  relationship. 


f^FDTD- 


J_,  f^l  2  ]  ^  3.75  X  IQ-^Sf  2  ] 
AzV  Ro  [l-Rl  Az  ll-pj 


(4) 


The  reflectance  and  conductivity  of  pure  gold  depend  on  the  morphology 
of  the  deposited  film.  Use  of  a  published  bulk  conductivity  [17]  of  2.44 
gfl-cm  in  equation  (3),  at  a  frequency  corresponding  to  a  wavelength  of 
10.0  pm,  yields  a  reflectivity  of  98.11  percent  and  a  synthetic  conductivity 
of  1.32  X  10^  S/ m  for  Az  =  0.3  pm.  On  the  other  hand,  published  optical 
properties  [11]  yield  a  reflectance  of  99.02  percent  and  real  and  synthetic 
conductivities  of  1.40  x  10^  S/m  and  2.56  x  10^  S/m,  respectively.  Because 
these  represent  pure  gold,  they  are  an  upper  limit  to  the  conductivity  and 
reflectance  that  can  be  achieved  in  a  real  QWIP.  The  actual  conductivity 
might  be  significantly  lower,  especially  if  alloyed  contacts  are  used. 
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Other  material  properties  used,  shown  in  table  1,  were  based  on  pub¬ 
lished  values  as  well  as  fits  to  experimental  data.  The  QW  stack  is  com¬ 
posed  of  many  alternating  layers  of,  for  example,  GaAs  and  AlGaAs; 
however,  since  the  layer  thicknesses  are  much  less  than  a  wavelength,  the 
stack  was  modeled  as  a  homogeneous  material  with  spatially  averaged 
properties.  The  permittivity  of  the  QW  stack  was  deduced  from  pub¬ 
lished  data  for  GaAs  and  AlAs  and  from  the  peak  wavelength  of  earlier 
grating-coupled  QWIPs.  The  conductivity  perpendicular  to  the  QW  stack 
was  deduced  from  Brewster-angle  absorption  measurements  on  unproc¬ 
essed  MBE  wafers.  Although  the  frequency-dependent  conductivity  of 
the  QW  stack  can  be  incorporated  into  the  FDTD  calculations,  we  used  a 
constant  conductivity  for  all  results  presented  here  so  that  the  coupling 
efficiency  could  be  seen  over  a  broad  spectral  range. 

Two  types  of  FDTD  runs  were  used.  The  most  straightforward  were  the 
steady-state  runs,  in  which  the  incident  field  on  the  pixel  was  a  sinusoid 
at  a  single  frequency.  In  this  case  the  FDTD  program  was  allowed  to  run 
until  the  fields  reached  a  steady  state.  The  XFDTD  program  was  set  up  to 
output  the  absorption  at  each  point  in  the  QW  stack  so  that  the  total  ab¬ 
sorption,  and  hence  the  QF,  could  be  computed.  Although  simple,  this 
technique  yields  the  QF  at  only  a  single  wavelength. 

In  the  second  technique,  the  incident  field  on  the  pixel  was  a  Gaussian- 
modulated  sinusoid,  as  shown  in  figure  4a.  The  Gaussian  modulation  in 
the  time  domain  corresponds  to  a  Gaussian  envelope  of  incident  frequen¬ 
cies  in  the  frequency  domain.  Since  the  FDTD  method  (and  all  of  our  ma¬ 
terials)  are  linear,  the  different  frequencies  in  the  incident  field  propagate 
independently  through  the  time  stepping.  Therefore,  as  shown  in  figure  4, 
Fourier  analysis  on  the  time-dependent  absorption  at  each  point  in  the 
QW  stack  yields  the  absorption  and  QF  at  all  of  the  incident  frequencies. 
In  other  words,  we  get  a  full  spectrum  of  QF  in  roughly  the  same  time 
that  the  steady-state  method  yields  the  value  at  a  single  wavelength. 


Table  1.  Material 
properties  used  in 
FDTD  calculations. 


Parallel  to  QWs 

Perpendicular  to  QWs 

Material 

Permittivity 
(Siemens /m) 

Gonductivity 

Permittivity 

(Siemens/m) 

Gonductivity 

GaAs  contact 

9.92 

0.0 

9.92 

0.0 

QW  stack 

9.92 

0.0 

9.92 

462 

Gold 

1.0 

1.0  X  10^  (“^l 

1.0 

1.0  X  10^  (“^l 

ZnS 

5.52 

0.0 

5.52 

0.0 

Epoxy 

2.25 

0.0 

2.25 

0.0 

(*)  "Synthetic"  conductivity  (see  text). 
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Figure  4.  Time-dependent  electric  fields  (left  column)  and  incident  and  absorbed  power  and 
QE  (right  column)  near  center  of  50-|Xm  pixel  with  3-|Xm  period  linear  grating  along  j/-direction: 
(a)  incident  field  (Gaussian-modulated  35.3-THz  sinusoid),  (b)  induced  for  j/-polarized 
incident  field,  and  (c)  induced  for  x-polarized  incident  field. 


The  main  problems  with  the  Gaussian-modulated  approach  are  practical 
ones.  Saving  the  full  set  of  time-dependent  fields  at  all  grid  points  in  the 
QW  stack  generates  huge  amounts  of  data.  For  example,  using  the 
XFDTD  "field  snapshot"  files  to  output  the  fields  at  all  time  steps  gener¬ 
ates  over  22  GB  of  data.  Writing  the  output  files  also  slows  the  program 
considerably.  However,  it  is  not  necessary  to  save  data  at  all  time  steps. 
First,  the  time  step  that  is  required  for  stability  in  the  FDTD  method  is 
highly  oversampled  relative  to  that  required  by  the  Nyquist  criterion.  For 
example,  the  Courant  timestep  required  for  stability  at  frequency/in  a 
material  with  index  n  and  with  cell  size  equal  to  1/10  the  wavelength  is 


At  = 


1 

lO/h'/S  ' 


(5) 
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which  is  5n-/3  times  smaller  than  the  Nyquist  sampling  interval  of  1  /  (2/) 
for  that  frequency.  By  sampling  the  fields  at  the  Nyquist  frequency  we 
reduce  the  amount  of  data  by  a  factor  of  ~27  for  GaAs  with  n  =  3.15. 
Second,  the  Gaussian  modulation  is  usually  set  up  with  a  relatively 
narrow  band  of  frequencies  in  the  incident  field  (corresponding,  for 
example,  to  a  wavelength  range  of  7-10  pm).  The  narrow  bandwidth  can 
be  exploited  by  further  undersamplmg  (below  the  Nyquist  frequency) 
such  that  the  frequency  band  is  aliased  down  to  lower  frequencies.  If  the 
undersampling  ratio  is  chosen  correctly,  the  shape  and  amplitude  of  the 
original  spectrum  can  be  preserved  after  the  aliasing.  By  combining  both 
techniques,  we  can  undersample  by  a  factor  as  high  as  -130  to  140  and 
still  recover  the  QE  at  all  of  the  frequencies  in  the  incident  field. 


Results 


Two-Dimensional,  Square-Cavity  Grating 

To  facilitate  comparison  with  earlier  modeling  and  with  experimental 
results,  we  computed  the  spectral  QE  for  a  conventional  QWIP  with  a 
3-pm-period  square  grating  using  the  EDTD  method.  As  shown  in  figure 
5,  the  structure  was  an  isolated  39-pm  pixel  with  40  QWs  surrounded  by 
9  pm  of  "epoxy."  The  grating  used  2.1-pm  square  cavities  that  were 
0.6  pm  deep.  A  cell  size  of  0.3  x  0.3  x  0.3  pm  was  used  with  8064  time 
steps.  The  time  step  length  was  equal  to  the  Courant  value  of 
5.78  X  10“^^  s.  The  electric  fields  had  decayed  to  less  than  1  percent  of 
their  peak  values  by  the  end  of  the  time  stepping. 

Eigure  6  shows  the  resultant  QE  of  the  structure  for  different  values  of  the 
grating  metal  EDTD  conductivity.  Several  features  are  noteworthy.  Eirst, 
the  high-conductivity  (s  =  1  x  10^  S/m)  curve  is  close  to  the  curve  for  PEC 
grating  metal,  verifying  that  the  synthetic  conductivity  method  gives  the 
appropriate  behavior  at  the  high-conductivity  limit. 

Second,  the  portion  of  the  high-  and  perfect-conductor  curves  above  8  pm 
are  very  similar  to  the  curves  shown  in  figure  7  of  Andersson  and 
Lundqvist  [4],  based  on  the  MEM  with  PEC  grating  metal.  The  similarity 
includes  the  double-humped  main  peak  between  9  and  10  pm,  as  well  as 
the  smaller  peak  between  8.5  and  9.0  pm.  This  similarity  exists  even 
though  the  layer  thicknesses  and  cavity  sizes  were  not  identical. 


Figure  5.  Cross  section 
of  square-grating- 
coupled  structure. 
Grating  is  3-pm  period 
with  2.1  X  0.6  pm 
cavities. 
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Figure  6.  Calculated 
spectral  QE  for  39-pm 
pixel  with  3-pm- 
period,  two- 
dimensional,  square- 
cavity  grating. 
"Synthetic" 
conductivity 
corresponding  to 
indicated  reflectances 
is  shown.  QW 
absorption  strength  is 
constant. 


-CT  =  1  X  10^  s/m  (R  =  75.0%) 
“a  =  3x10*  S/m  (R  =  91.7%) 
-a  =  1  X  10*  S/m  (R  =  97.5%) 
-a  =  1  X  10'  S/m  (R  =  99.8%) 
-a  =  PEC  (R  =  100%) 


Figure  7.  Operation  of  spiral  antenna:  (a)  operating  regions  and  (b)  field  distributions  in  transmission 
line  region  with  and  without  etching. 


Third,  the  QE  is  significantly  degraded  by  finite  conductivity,  even  for 
conductivities  corresponding  to  R  =  97.5  percent.  Note  particularly  the 
suppression  of  the  long-wavelength  peak  at  about  9.65  pm.  This  peak  is 
beyond  the  9.45-pm  resonance  of  the  grating  and  therefore  corresponds  to 
evanescent  modes.  The  evanescent  modes  (and  near-resonance  propagat¬ 
ing  modes)  that  couple  to  the  QW  stack  have  electric  field  perpendicular 
to  the  grating  surface  and  magnetic  field  along  the  surface.  The  tangential 
magnetic  fields  imply  surface  currents  in  the  grating  metal  that  are  sup¬ 
pressed,  along  with  the  evanescent  modes  themselves,  when  the  grating 
conductivity  is  reduced.  We  note  that  the  curve  for  R  =  91.7  percent  ap¬ 
proximately  matches  the  shape  and  magnitude  of  many  real  QWIP  detec¬ 
tors,  suggesting  that  grating  metal  conductivity  may  be  a  major  factor 
that  limits  the  performance  of  these  type  of  QWlPs. 
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Spiral  Antenna 

The  spiral  antenna  offers  a  number  of  attractive  features  for  use  as  an  op¬ 
tical  coupler  in  QWIPs.  It  is  well  known  for  wide  spectral  bandwidth, 
typically  8:1  or  even  32:1  in  the  microwave  region  [18].  This  is  much 
wider  than  the  2:1  or  3:1  ratio  required  for  multicolor  QWIPs.  (However, 
use  of  a  backplane  reflector  reduces  the  bandwidth.)  Finally,  the  field  pat¬ 
terns  around  a  two-  or  four-arm  spiral  antenna  have  large  components 
perpendicular  to  the  antenna  plane,  which  should  strongly  couple  to  a 
closely  placed  QW  stack. 

Figure  7a  shows  the  basic  operating  regions  of  a  standard  spiral  antenna. 
In  the  so-called  radiation  zones,  the  currents  in  adjacent  windings  are  in 
phase  and  lead  to  radiation  of  power  from  the  antenna  (or  absorption  of 
power  from  an  incident  field).  In  the  transmission-line  regions,  the  cur¬ 
rent  in  adjacent  windings  is  approximately  180°  out  of  phase  so  that 
power  is  carried  along  the  windings  as  in  a  standard  two-wire  transmis¬ 
sion  line.  As  indicated  in  figure  7b,  the  electric  fields  in  the  transmission¬ 
line  regions  have  large  components  perpendicular  to  the  plane  of  the  an¬ 
tenna  that  can  couple  to  a  closely  spaced  QW  stack.  The  QW  stack  can  be 
thought  of  as  loading  the  antenna  so  that  the  reradiated  power  is  less  than 
the  incident  power.  This  power  "loss"  from  an  antenna  perspective  is  the 
desired  detection  from  a  QWIP  perspective. 

The  coupling  strength  between  the  antenna  modes  and  the  QW  stack  can 
be  further  increased  if  the  QW  material  between  the  windings  is  removed 
and  replaced  with  a  material  with  a  lower  refractive  index  (such  as 
vacuum  or  ZnS).  As  shown  in  figure  7b,  the  electric  field  lines  tend  to 
flow  through  the  higher  index  QW  material,  thereby  extending  the  cou¬ 
pling  farther  into  the  QW  stack.  This  removal  of  about  half  of  the  QW 
material  also  lowers  the  dark  current  by  the  same  ratio. 

Figure  3  shows  the  antenna  pattern  that  was  used  for  this  work.  The  unit 
cell  is  a  single,  four-arm  spiral  that  is  replicated  in  a  two-dimensional  ar¬ 
ray  to  fill  the  pixel.  The  individual  spirals  are  interconnected  so  that  a 
single  via  hole  through  the  backplane  dielectric  can  be  used  to  make  con¬ 
tact  to  the  entire  array  of  antennas,  which  also  acts  as  one  of  the  electrical 
contacts  to  the  detector.  The  square  spiral  was  used  to  ease  fabrication  us¬ 
ing  electron-beam  lithography  (EBL).  The  lines  and  spaces  are  about 
0.32  pm  wide,  which  is  easily  achievable  with  EBL.  A  cell  size  of  0.16  x 
0.16  X  0.16  pm  was  used  with  19,328  time  steps.  The  time  step  length  was 
equal  to  0.75  times  the  Courant  value,  or  2.41  x  lO'^*’  s.  The  electric  fields 
had  decayed  to  less  than  1  percent  of  their  peak  values  by  the  end  of  the 
time  stepping.  We  modeled  a  3  x  3  array  of  spirals  that  includes  the  inter¬ 
actions  between  neighboring  spirals  and  would  be  appropriate  for  a 
20-pm  pixel. 
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Figure  8  shows  the  calculated  spectral  QE  from  the  3x3  array  for  struc¬ 
tures  with  QW  stack  thicknesses  of  0.33  and  0.83  pm,  corresponding  to  5 
and  12  QWs,  respectively.  The  QEs  of  20  to  30  percent  for  12  QWs  and  10 
to  20  percent  for  5  QWs  are  relatively  high,  considering  the  small  num¬ 
bers  of  QWs.  Eurthermore,  since  the  dark-current-limited  D*  is  propor¬ 
tional  to  QE  '  where  is  the  number  of  wells  and  is  the 

dark  current,  a  QE  of  20  percent  with  12  QWs  (or  13  percent  with  5  QWs) 
and  a  50-percent  dark  current  reduction  will  yield  the  same  D*  as  52 
percent  QE  with  40  QWs  and  full  dark  current.  The  figure  also  shows 
curves  for  a  synthetic  conductivity  of  10^  S/ m.  The  QE  is  substantially 
reduced  below  the  PEC  curves,  but  as  previously  stated,  the  physical 
conductivity  implied  by  this  synthetic  conductivity  is  not  clear  for  a  layer 
only  one  cell  thick. 

Eigure  9  shows  the  absorption  patterns  in  the  QW  stack  near  the  spiral  for 
two  different  wavelengths  of  illumination.  The  pattern  oscillates  as  one 
moves  along  the  winding,  as  expected  from  the  previously  stated  operat¬ 
ing  mechanism.  The  physical  mechanism  behind  the  antenna's  broad 
spectral  bandwidth  is  also  apparent.  Note  that  the  absorption  pattern  ro¬ 
tates  around  the  antenna  windings  as  the  wavelength  is  varied,  but  main¬ 
tains  nearly  the  same  total  absorption. 


Figure  8.  Calculated 
spectral  QE  for 
antenna-coupled 
QWIP  with  3x3  array 
of  spiral  antennas. 
Conductivity  is 
"synthetic" 
conductivity 
corresponding  to 
indicated  reflectances. 
QW  absorption 
strength  is  constant. 


12  QW,  3x3,  PEC 
5QW,  3x3,  PEC 


5  QW,  3  x3,  o  =  10®  S/m 
12  QW,  3x3,  o  =  10®  S/m 


Wavelength  (|im) 


Figure  9.  Absorption 
in  QW  stack  near 
spiral  antennas  at 

(a)  9.42  pm  and 

(b)  9.06  pm.  Structure 
has  12  QWs  with  PEC 
antenna. 
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Jigsaw  Antenna 

The  jigsaw  antenna  pattern  is  shown  in  figure  10.  It  was  originally  de¬ 
signed  as  a  set  of  vertical  and  horizontal  lines  that  had  U-shaped  shunts 
to  induce  out-of-phase  currents  in  neighboring  lines  similar  to  those  in¬ 
duced  in  the  spiral  antenna,  but  with  a  longer  decay  length  away  from 
the  antenna  plane.  The  pattern  could  then  couple  to  thicker  QW  stacks  for 
applications  where  maximum  QE  is  the  most  important  criterion. 

However,  after  optimizing  the  dimensions  using  the  infinite-pixel  model, 
we  found  that  the  pattern  contains  an  overall  period  of  ~18  pm  and  a 
"sub-period"  (of  the  lines  in  the  pattern)  of  ~3  pm.  Therefore,  the  pattern 
acts  much  like  a  3-pm  period  planar  metal  grating.  However,  since  the 
pattern  has  an  18-pm  period,  it  supports  approximately  36  (18/3, 
squared)  times  as  many  propagating  modes  as  the  pure  3-pm  period  grat¬ 
ing.  With  the  proper  layer  thicknesses,  these  modes  combine  to  suppress 
rediffraction  out  of  the  detector,  thereby  yielding  high  QE. 

We  modeled  a  jigsaw  antenna-coupled  pixel  with  40  QWs  that  contained 
2x2  unit  cells,  3x3  unit  cells,  and  an  infinitely  periodic  pixel.  A  cell  size 
of  0.3  X  0.3  X  0.3  pm  was  used  with  16,384  time  steps.  The  time  step  length 
was  equal  to  the  Courant  value  of  5.65  x  10“^^  s.  The  results  are  shown  in 
figure  11.  The  figure  includes  results  for  both  PEC  and  finite-conductivity 
antennas,  but  with  the  previously  mentioned  warning  about  the  uncer¬ 
tainty  of  the  synthetic  conductivity  method  for  a  one-cell-thick  antenna. 
Note  that  the  spectral  QE  is  high  over  a  broad  spectral  range  for  the 
infinite  pixel.  However,  the  QE  of  the  3x3  and  2x2  period  pixels  are 
progressively  degraded  due  to  finite-pixel  size  effects,  consistent  with  this 
design's  use  of  propagating  modes. 


Figure  10.  2  x  2 
periods  of  jigsaw 
anteirna  in  38-pm 
pixel.  Metal  is 
dark. 


38-|im  pixel 


18-(im  unit  cell 
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Figure  11.  Calculated 
QE  for  several  pixel 
sizes  using  jigsaw 
antenna  coupler. 
Finite  pixels  are 
calculated  by  FDTD. 
Infinite  pixel  is 
calculated  by  MEM. 
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Conclusions 


Planar  antennas  offer  several  benefits  as  optical  couplers  in  QWIPs.  Use 
of  planar  metal  layers  eliminates  the  need  for  a  grating  etch  and  the  po¬ 
tential  nonuniformity  associated  with  that  etch.  All  of  the  processing  ex¬ 
cept  for  deposition  of  the  AR  coating  is  done  on  the  front  (MBE  growth) 
side  of  the  wafer  so  there  is  no  need  for  backside  processing  of  individual 
die  or  handling  of  fully  thinned  wafers,  as  with  some  other  approaches. 

As  demonstrated  by  the  calculated  spectral  QE  for  the  spiral  antenna  ar¬ 
ray,  a  relatively  high  QE  can  be  obtained  over  a  wide  spectral  band,  even 
with  a  reduced  number  of  QWs.  This  should  permit  increased  detectivity 
in  dark-current-limited  QWIPs  at  elevated  operating  temperature. 

We  have  also  presented  the  first  numerical  modeling  of  QWIPs  using  the 
PDTD  method  and  demonstrated  the  utility  of  the  technique  for  predict¬ 
ing  QE  for  finite-sized  pixels.  The  undersampled  spectral  QE  technique  is 
particularly  useful,  since  it  yields  the  full  spectral  response  in  a  single 
EDTD  run.  After  more  validation  of  the  predicted  results  with  experi¬ 
ment,  the  technique  should  be  a  very  useful  tool  for  designing  and  opti¬ 
mizing  realistically  sized  QWIP  pixels. 
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